Title : will be set by the publisher 
Editors : will be set by the publisher 
EAS Publications Series, Vol. ?, 2008 



MULTIMODE RADIATIVE BETA CEPHEI AND DELTA 
CEPHEI MODELS 

Smolec R. 1 



Abstract. As a by-product of f3 Cephei and 8 Cephei radiative model 
surveys we have found interesting cases of multimode pulsation. In 
case of p Cephei stars, we found two multimode domains with two or 
three modes being involved. The origin of the multimode pulsation can 
be traced to one of the two different mechanisms: either to the non- 
resonant coupling of the two excited modes (double-mode pulsation) 
or to the 2wi ~ uoo + u>2 parametric resonance (triple-mode pulsation) . 
In case of S Cephei models, a triple-mode pulsation domain, connected 
with the 2:1 resonance, 2ll>o ~ u>2, is found. These models are of theo- 
retical interest only, since they do not obey a proper mass-luminosity 
relation. 



1 Introduction 

Doubly-periodic classical pulsators (beat Cepheids, beat RR Lyrae stars) are the 
simplest and one of the most important examples of multimode stellar pulsation. 
Observed stellar variation results from simultaneous oscillation in two low-order 
radial modes of pulsation. More than 30 years ago it was realized, that simple 
linear theory may be used to calculate the masses of such pulsators (Petersen 
I1973p . However it took many years to understand and to model doubly-periodic 
pulsators. Some questions remain open even today. 

Throughout this paper we will focus on doubly-periodic pulsations. We stress 
that doubly-periodic is not equivalent to double-mode, as more than two pulsation 
modes may be involved in doubly-periodic oscillations if resonant mode interaction 
is present. Such examples will be shown in this paper. 

Attempts to model the doubly-periodic pulsators started with early radiative 
lagrangian hydrocodes. These first nonlinear calculations failed to reproduce the 
doubly-periodic behaviour, but they were very incomplete, as direct nonlinear 
calculations are very time-consuming. Introduction and implementation of relax- 
ation technique by Baker & von Sengbusch ([1969]) and Stellingwerf (1974) allowed 
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for efficient search for multimode solutions. Relaxation technique allows for fast 
convergence to limit cycle (periodic, finite-amplitude pulsation) and provides in- 
formation on its stability through the Floquet stability coefficients (or so-called 
switching rates). Two Floquet coefficients are of interest, when one searches for 
doubly-periodic models with fundamental and first overtone modes simultaneously 
excited. 70,1 measures the stability of the first overtone limit cycle (switching rate 
toward fundamental mode), while 71.0 measures the stability of the fundamental 
mode (switching rate toward first overtone) . If Floquet coefficients of both limit 
cycles are positive, neither of the limit cycles is stable and simultaneous pulsation 
in both modes is unavoidable. 

Using relaxation technique Kovacs & Buchler (1988) performed detailed survey 
in search for the multimode pulsation in RR Lyrae models. They found that 
the models close to the 2luq — ^3 resonance center exhibit the doubly-periodic 
behaviour, however with periods inconsistent with observed ones. Physical reasons 
for resonant, doubly-periodic pulsations were analysed in detail by Dziembowski & 
Kovacs (1984) through the analysis of the amplitude equations. Linearly damped 
third overtone acts as an energy sink for the linearly unstable fundamental mode. 
As a result amplitude of the fundamental mode is lowered and otherwise dominant 
fundamental mode is no longer able to saturate the pulsational instability. This 
allows the growth of the first overtone mode. Resulting triple-mode solution is 
apparently doubly-periodic due to nonlinear frequency-lock phenomenon. In the 
nonlinear, resonant regime 0)3 is exactly equal to 2Qq and third overtone is not 
visible as a separate frequency, but contributes to the harmonic frequencies of 
the fundamental mode. Other examples of resonant triple-mode, doubly-periodic 
solutions will be presented in the next Sections. 

Kovacs & Buchler (1988) also found non-resonant double-mode models through 
playing with artificial viscosity parameters. Artificial viscosity provides unphysical 
dissipation necessary to handle the shocks in the radiative codes. Model proper- 
ties, specially the amplitudes depend on artificial viscosity (eg. Kovacs H990p . 
With artificial viscosity parameters chosen to satisfy the observed amplitudes no 
multimode solution was found. However models that are less dissipative due to 
reduced artificial viscosity, exhibit steady doubly-periodic, double-mode pulsation 
as was investigated in detail by Kovacs & Buchler (|1993p . Their models however, 
strongly depend on numerical details and do not reproduce all the observational 
constraints. Nevertheless, experiments with artificial viscosity motivated the de- 
velopment of numerically less dissipative codes with physical dissipation instead 
of artificial one. The success came with the introduction of turbulent convec- 
tion energy transfer recipies into the hydrocodes. Kollath et al. (|1998p reported 
double-mode Cepheid models, while double-mode RR Lyrae models were reported 
by Feuchtinger (|1998j) . These double- mode models originate from non-resonant 
coupling of the pulsation modes and reproduce the observed periods and period 
ratios. 

In this paper we present new interesting examples of purely radiative doubly- 
periodic models. In Section 2 we discuss numerically robust, resonant and non- 
resonant (3 Cephei models, while in Section 3 resonant 5 Cephei models will be 
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presented. 

2 Multimode radiative 3 Cephei models 

P Cephei stars are main sequence, early B-type pulsators. Acoustic modes and 
low-order gravity modes are driven through the K-mechanism acting in the iron 
opacity bump (e.g. Pamyatnykh 1999). Most of these stars are multiperiodic. 
Non-radial modes dominate, although radial pulsation is observed in several stars. 
Under some simplified assumptions Smolec & Moskalik (|2007|) studied the ampli- 
tude saturation in these stars, using purely radial, radiative codes. In case of these 
stars convection may be safely disregarded, as /3 Cephei stars are too hot to have 
an effective energy transfer by convection ( P amvat nvkh 1 1 9 9"9l Smolec & Moskalik 
I2007|) . Model properties are also independent of artificial viscosity, contrary to the 
classical pulsators case (see e.g. Kovacs [T990l and Smolec & Moskalik l2007p . Thus, 
calculated radiative models are very robust. As a by-product of this analysis sev- 
eral multimode models were found. The modal selection depends on the opacities 
being used. Here we discuss the models computed with the OPAL opacities, as 
they represent more interesting and clear case. Multimode models calculated with 
the OP opacities are discussed in detail by Smolec & Moskalik (2007). 

Computation of the Floquet stability coefficients for the models located on 
10M Q evolutionary track revealed two multimode domains where first overtone 
and fundamental mode are simultaneously unstable (both Floquet coefficients, 
7o,i and 71, o> are positive). Linear frequencies of the models show, that one of 
the domains is clearly connected with the high-order parametric resonance 2wi ~ 
ujq + cl>2- Situation is illustrated in Fig. [T] where Floquet coefficients and linear 
growth rates of the modes, 70 and 71, are plotted versus the A parameter, that 
characterizes the proximity of the resonance, A = 2lui/(luq + 102). Starting from 
the blue edge, switching rate toward fundamental mode is negative and thus, 
first overtone is stable. Then, exactly at the resonance center, a prominent peak 
of the switching rate toward fundamental mode is visible. Otherwise negative, 
70,1 becomes positive at the resonance center. Thus, the resonance leads to the 
destabilization of the first overtone limit cycle. Resulting triple-mode domain 
appears as an island in the middle of the first overtone pulsation domain. The 
resonance does not affect the stability of the fundamental mode, as is well visible 
in Fig. [TJ and as was shown by Kovacs & Buchler (|1993[) through the analysis of 
the amplitude equations. Second multimode domain located at A « I. Of (Fig. Q} 
is not connected with any resonance. Thus, only fundamental and first overtone 
modes are involved in double-mode oscillations. This domain separates the first 
overtone and fundamental mode pulsation domains as is expected in case of non- 
resonant mode coupling. 

We have also performed some nonlinear calculations to approach the doubly- 
periodic solution. Due to the very low growth rates, several hundred thousand 
cycles are necessary in nonlinear integrations to approach the attractor. Thus, 
we focused on only one model located in the non-resonant double-mode domain. 
Nonlinear integrations were initiated with four different initial conditions (velocity 
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Fig. 1. Linear growth rates (thin lines) and Floquet coefficients (thick lines) plotted 
versus the A parameter (A = 2wi/(wo + ^2)) for the sequence of (5 Cephei models. 
Shaded bar indicates the modal selection as derived from Floquet coefficients (IO - 
first overtone, F - fundamental mode, TM - triple-mode, DM - double-mode pulsation 
domains) . 



perturbations). Relative surface displacement, 8R/R, was then followed and the 
amplitudes of the fundamental and first overtone modes were extracted through 
the use of the analytic signal method (Kollath & Buchler 2001). Fig. [2] presents 
the results in the modal amplitude phase-space. Transient evolution is extremely 
slow and none of the trajectories reached the double-mode attractor. However its 
location is clearly visible: two trajectories evolve toward the attractor from the 
right, and two of them from the left. 

Both of the discussed multimode domains represent clear examples of radial 
multimodc pulsation. Origin of the multimode state has clear physical interpre- 
tation. Computed models are robust, independent of numerical details. Also 
radiative codes are fully applicable in case of (3 Cephei stars. Nevertheless, two 
radial modes simultaneously excited, were not observed in these stars up to date. 
These results however, shed some light on the reasons of multiperiodic nature of 
f3 Cephei stars. Saturation properties of the non-radial acoustic modes should be 
very similar to the properties of the radial modes. Thus, non-resonant coupling 
of non-radial modes may lead to multimode pulsations. Taking into account the 
dense spectrum of the non-radial modes, resonant coupling leading to multimode 
pulsation is also very likely (Dziembowski [20071 Smolec & Moskalik l2007p . 
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Fig. 2. Temporal evolution of the first overtone, Ai, and the fundamental mode, Ao, 
amplitudes of the non-resonant /3 Cephei model. Integrations were initialized with four 
different initial conditions and followed for more than hundred thousand cycles. 



3 Triple-mode radiative 5 Cephei models 

Triple-mode models we will discuss in this Section, are connected with the 2:1 res- 
onance between the linearly unstable fundamental mode and the linearly damped 
second overtone mode, 2ujq — L02- This resonance plays a crucial role in shap- 
ing the Hcrtzsprung bump progression in fundamental mode Cepheids (Kovacs & 
Buchler [H)g9l Buchler, Moskalik & Kovacs [P9"9"0)) . These so-called bump Cepheids 
are indeed a double-mode pulsators. However, due to frequency-lock phenomenon, 
second overtone is not visible as a separate frequency, but manifests as a nonlinear 
distortion of the light /velocity curves. Observationally, resonance center falls at 
fundamental mode periods around Pq 10 days. Due to the resonant interaction 
with the damped second overtone, amplitude of the fundamental mode is lowered, 
which may lead to the destabilization of its limit cycle. Indeed, Buchler et al. 
(1997) show that for some models, 7^0 > 0, close to the resonance center and 
thus, fundamental mode becomes unstable. However, these models are located 
relatively close to the blue edge and the first overtone is stable (70,1 < 0). Thus, 
the star pulsates in the first overtone, and simultaneous pulsation in fundamental 
and first overtone mode is not possible. This is in agreement with observations, 
since double- mode first overtone/fundamental mode Cepheids have much shorter 
periods and, as linear calculations show, lie in a non-resonant regime. 

Using different mass-luminosity relation we have shifted the resonance center 
toward much lower temperatures and longer fundamental mode periods. As a re- 
sult resonance center falls at the temperature range, where first overtone is already 
unstable, 70,1 > 0. Situation is illustrated in the left panel of Fig. (3) For higher 
temperatures both limit cycles are stable. Thus, pulsation in either fundamen- 
tal mode or first overtone is possible. For lower temperatures first overtone limit 
cycle becomes unstable and fundamental mode pulsation domain begins. Then, 
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close to the resonance center, otherwise stable, fundamental mode becomes un- 
stable. This is a direct consequence of lower amplitude of fundamental mode, as 
is shown in right panel of Fig. [3J Amplitude is lowered because part of the fun- 
damental mode energy is resonantly transferred to the damped second overtone. 
Triple-mode pulsation, with two frequencies (second overtone is phase-locked) is 
possible. We stress that these models are of theoretical interest only, as they result 
from arbitrary mass-luminosity relation, that is inconsistent with observations. 




Fig. 3. Floquct coefficients plotted versus the model's temperature (left panel). Shaded 
bar indicates modal selection derived from Floquet coefficients and from nonlinear in- 
tegrations. Destabilization of the fundamental mode limit cycle (71,0 > 0) is clearly 
connected with its lower amplitude (right panel). Arrows indicate the position of the 
2cjo — 0J2 resonace center. 



Several nonlinear integrations were performed to approach the triple-mode at- 
tractor. Relative surface displacement of the model, SR/R, was followed. Analytic 
signal method was used to extract the amplitudes, Sq and SI, corresponding to the 
variability with frequencies loq and ui\. S\ is directly related to the amplitude of 
the first overtone pulsation, S\ = A\. However, due to resonant coupling between 
the fundamental and second overtone modes, Sq contains contribution from both 
of the modes, namely: 

So = 3?(ao + ^21^0^2) 

where ao and 02 are complex amplitudes corresponding to fundamental and second 
overtone modes respectively and ft.21 is complex coefficient (Kovacs & Buchler 
1989). We note that the topology of the S0-S1 phase-space is the same as the 
topology of the Aq-A% phase-space. Exemplary trajectories in S0-S1 phase-space 
are shown if Figs. [4] and [5] In Fig. [4] we show the trajectories for five models with 
different temperatures. Each model on this Figure was kicked with the same initial 
condition. We also investigated the details of the phase-space by integrating the 
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Fig. 4. Phase-space plots for models with different temperatures. All integrations were 
started with the same initial condition and continued till the attractor was reached. 

same static model with different initial conditions. Results are shown in Fig. [5] for 
two different models: 4690K model (left panel) and 4720K (right panel). Several 
interesting features are visible in Figs. [4] and O 

(i) The 4730K model trajectory evolves toward single-frequency solution. At 
lower temperature (4720K) triple-mode attractor appears. Trajectories for the 
4690K model (Figs. 2] and \5\ left) evolve toward a pronounced spiral attractor. 
Such attractor is not possible in non-resonant case. For lower temperatures spiral 
becomes less pronounced and disappears for 4660K model. For slightly cooler 
model triple-mode attractor also disappears. 

(ii) For lower temperatures triple-mode solution is the only attractor of the 
system (Fig. [S] left). For these models both Floquet coefficients are positive (Fig. [3] 
left). 

(hi) For higher temperatures (Fig. [5] right) hysteresis is possible. Depending 
on the initial conditions, trajectories evolve toward two different attractors: a 
triple-mode attractor and a single-frequency fundamental mode attractor. For 
these models fundamental mode is stable (71,0 < 0, Fig. [3]). Similar hysteresis 
(double- mode/fundamental mode) is present in non-resonant convective Cepheid 
and RR Lyrae models (Kollath & Buchler 2001). However here, the single- 
frequency attractor is a double-mode one, being a bump Cepheid. 
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Fig. 5. Structure of the phase-space for 4690K model (left panel) and 4720K model 
(right panel). In the latter case two attractors are visible. 
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